Possible Coexistence of Rotational and Ferroelectric Lattice Distortions in 
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The competitions between ferroelectric and rotational in- 
stabilities in rhombohedral PhZr^Tii-^Oa near x = 0.5 are 
investigated using first principles density functional supercell 
calculations. As expected, we find a strong ferroelectric insta- 
bility. However, we also find a substantial R-point rotational 
instability, close to but not as deep as the ferroelectric one. 
This is similar to the situation in pure PbZrC>3. These two 
instabilities are both strongly pressure dependent, but in op- 
posite directions so that lattice compression of less than 1% 
is sufficient to change their ordering. Because of this, local 
stress fields due to B-site cation disorder may lead to coex- 
istence of both types of instability are likely present in the 
alloy near the morphotropic phase boundary. 

Rhombohedral PbZr x Tii_2;03 (PZT) ceramic alloys 
with compositions near the morphotropic phase bound- 
ary (MPB) around x — 0.52. form the basis of most piezo- 
electric transducer devices .El This is due to a combination 
of high response, realizable strains in the tenths of 
range and favorable weak temperature dependencies. 
Recently, there has been renewed scientific interest in 
these materials, driven partly by the discovery of nearly 
order of magnitude higher performance in related relaxor 
single crystals,Q'd and partly by the fundamental under- 
standing of them, being obtained by experimental inves- 
tigation and first principles calculations. 

The microscopic physics underlying the low tempera- 
ture phases of the end-points is well known. In PbTiC>3 
ferroelectricity is due to condensation of a T\s unsta- 
ble phonon where the oxygen octahedra shift against the 
cations. The ground state structure has shifts along (001) 
with a tetragonal lattice strain that stabilizes this direc- 
tion. A rhombohedral ferroelectric (FE) phase with (111) 
shifts is not favored and does not occur because of the 
large electronic hybridization between Pb and O, as may 
be seen by comparison with .BaTiOa, which has a rhom- 
bohedral FE ground state-ETO 

PhZcD 3 has a complex anti-ferroelectric ground 
statetrO that may be viewed as arising from the ideal 
cubic perovskite structure bynra, combination of strong 
zone boundary instabilities.Ejii-3 Significantly, the fer- 
rodistortive mode that gives the FE ground state 
of PbTiC>3 is also found strongly unstable in density 
functional (DF) studies of cubic perovskite PbZrC>3 even 
though the ground state is not FE. Also the instabili- 
ties of cubic perovskite PbZrC>3 are much stronger than 



for P_UTip3, with energies of 0.20 - 0.25 eV per formula 
unit]i3'tLa The actual structure arises from a delicate bal- 
ance between modes characterized as octahedral rota- 
tions, octahedral distortions and off-centering. With the 
addition of small amounts of Ti, the PZT phase diagram 
shows a transition to FE behavior, in other words freezing 
in of the Tis instability. However, the actual situation is 
no doubt more complex than this. First of all, in the Zr 
rich part of the phase diagram, there is a low temperature 
rhombohedral phase in addition to the high temperature 
phase characterized by the pure Tis displacement. This 
low temperature phase is associated with a co-existence 
of the frozen in r^-EE-instability with rotations of the 
oxygen octahedraEOEZl Secondly, local probes indicate 
a complex local structure that differs substantially from 
the average diffraction structure in Zr rich PZT.atid--also 
in the related Pb based relaxor single crystals .B'ta~E3 

Complete phonon dispersions of cubic perovskite 
PbTiC>3 and PbZrC>3 determined by DF calculations were 
reported by Ghosez et al.E3 Both materials show both 
FE (ris) and rotational (R25 type) instabilities, though 
in the titanate the rotational (R25) instability is weak 
and occurs only in small regions of the zone around the 
R and M points. As mentioned, the R25 unstable mode 
consists of rotation of the oxygen octahedra around the 
transition metal. The rapid upward dispersion away from 
R and M is notable; it reflects the rigidity of the octahe- 
dra. In cubic perovskite structure PbZrOs the R25 mode 
is more unstable and disperses upwards from R weakly, 
implying more deformable octahedra, ^a. point that was 
associated with pressure dependencies.LJ It should, how- 
ever, be noted that in materials like PbZrC>3 where the 
ideal structure is highly unstable, the size of the various 
instabilities is not simply related to the magnitudes of 
the corresponding imaginary phonon frequencies in the 
cubic structure. For example, the phonon frequencies in- 
dicate that the R and M point instabilities associated 
with octahedral rotation in P. bZrQ 3 are much more un- 
stable than the FE Ti5 mod e. aH but, energetically, the 
FE mode is more unstable.HLj 

Here, we use a simple ordered supercell with composi- 
tion x — 0.5 to investigate the relative strengths of the 
FE and rotational instabilities in this region. We find 
that the FE instability is stronger, as expected, but only 
marginally so, and that the balance between these may 
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be easily switched by modest strains at the sub 1% level. 
Noting the difference in ionic radii of Ti and Zr (and the 
cell volumes of PbZrC>3 and PbTiOa) and that there is. 
no evidence for cation ordering in PZT near the MPB,E^I 
one may conclude, first of all, that these two types of 
instability may coexist in alloys near the MPB and, sec- 
ondly, that it may be helpful to build rotational degrees 
of freedom into effective Hamiltonians. 

We focus on the rhombohedral side of the MPB ne- 
glecting the rhombohedral strain that is kapwn to be very 
small in contrast to the tetragonal phaseE-3 The supercell 
used to model the alloy is a 10 atom FCC cell with alter- 
nating Ti and Zr layers along the [111] direction. This 
is .the same as one of the cells used by Saghi-Szabo et 
al.u to investigate piezoelectricity on the tetragonal side 
of the MPB-ajid. Ramer et al. to investigate microscopic 
stress nelds.c3't£l The present DF calculations were done 
within the Hcdin-Lundqvist local density approximation 
(LDA) using the general potential linearized augmented 
plane-wave methocO with local orbital extensionaJ.0 re- 
lax linearization errors and treat semi-core stateso The 
Brillouin zone samplings were done using 4x4x4 special 
k point meshes (note this is for the doubled perovskite 
cell). Well converged basis sets of over 1650 functions 
were used with sphere radii of 2.0, 1.83, 2.25 and 1.47 
ao for Zr, Ti, Pb and O, respectively. Forces were cal- 
culated by the method of Yu et al.E3 Calculations were 
done for three lattice parameters in 1% increments, i.e. 
7.555, 7.631 and 7.708 a . While the ideal low temper- 
ature cubic perovskite structure cannot be accessed ex- 
perimentally, Jaffe et alEl report an effective (by cell 
volume in the tetragonal phase) value of 7.73 oo at room 
temperature, while Noheda et alr3 obtain 7.69 ao for 
both the tetragonal (slightly above room temperature) 
and monoclinic (20 K) phases at a slightly more Zr rich 
composition, x — 0.52. Within the LDA we obtain a 
slightly smaller effective lattice parameter for the super- 
cell: 7.55 oo, with cubic symmetry (but including the O 
breathing), and 7.59 ao with a full relaxation into the 
ferroelectric structure. Such 1 — 2% smaller lattice pa- 
rameters relative to experiment are typical of LDA errors 
in this class of materials S3c3 

In the ideal cubic structure, the supercell has one inter- 
nal parameter corresponding to breathing of the octahe- 
dron around Ti. At a = 7.631 ao relaxation of this gives 
an energy gain of 33 mRy (all energies are for the 10 atom 
cell) with a 0.116 ao reduction of the Ti-0 bond lengths 
- very close to what would be expected (0.11 ao) based 
on the difference in Ti and Zr ionic radii. This isotropic 



TABLE I. Calculated supercell vibrational frequencies 
(cm -1 ) of phonons compatible with rhombohedral symmetry. 
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FIG. 1. Relaxation along the FE unstable mode (energies 
for 10 atom supercells) for the reference volume (circles), 1% 
compression (squares), and 1% expansion (diamonds). The 
instability decreases under pressure. 

breathing corresponds to the highest phonon branch as 
shown in Table |. The frequency of 805 cm -1 is higher 
than that calculated in either pure PbZ*03 or PbTiOs 
(725 cm -1 and 612 cm -1 , respectively)^! as can be ex- 
pected from bond length considerations (note that the 
O is between a Zr and Ti yielding a short relaxed Ti-0 
bond length). 

We computed the remaining phonons compatible with 
a rhombohedral symmetry (these are the Tis and zone 
folded R15 modes of the simple cubic perovskite) from 
atomic forces for a variety of small distortions about the 
breathed FCC structure. The dynamical matrix was de- 
termined by least squares fit to these, and then diagonal- 
ized to obtain the phonon eigenvectors and frequencies 
(given in Table |) for the three volumes. Additionally, 
we calculated the energetics of the simple perovskite R25 
rotational mode (which, however, is not strictly speak- 
ing a true pure mode for the FCC supercell). We note 
that the ferroelectric and rotational modes and the Ti5 
and folded R15 modes belong to different symmetries and 
therefore do not mix at the harmonic level {e.g. the fer- 
roelectric mode breaks inversion). 

The 500 cm -1 mode in Table | involves mainly O dis- 
placements, like the breating mode: small octahedra tilt- 
ing mixed with stretching along [001]. The other modes 
above 300 cm -1 involve distortions of the O octahedra, 
while the lowest stable mode involves mainly transition 
metal off-centering in the octahedra. In rhombohedral 
symmetry there is also one marginally unstable and one 
unstable mode - the R15 mode and the FE Tis insta- 
bility, respectively. In the ideal perovskite cell and the 
true disordered alloy, they belong to the same phonon 
branch and a rough interpolation suggests that the FE 
instability extends over the entire Brillouin zone as ex- 
pected by comparison with the full dispersion curve in 
PbZrOa (Ref. |2l]). Energy minimization along the FE 
eigenmode at a = 7.631 ao yields an energy gain of 11 
mRy (again per doubled cell). The local minimum com- 
patible with rhombohedral symmetry, however, occurs 
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away from this line and the relaxation provides futhcr 
energetic gain of 7 mRy. Its coordinated compare well 
with the displacement pattern obtained by Saghi-Szabo 
it et al. for the tetragonal FE state (note that at the 
harmonic level the tetragonal and rhombohedral eigen- 
vectors are the same).H As shown in Fig.jl], the instability 
is energetically disfavored by compression and favored by 
expansion. Over the range of 2% in lattice parameter we 
find a variation of the energy gain along the FE eigen- 
mode of about 5 mRy. Further, the variation in depth 
of the local minima are larger: 11, 18 and 22 mRy for 
a = 7.555 do, a — 7.631 clq and a — 7.708 ao respectively. 

The precise mechanism of the high electromechanical 
response in PZT and the related relaxor crystals is still 
not fully established but there are strong indications that 
polarization rotation from rhombohedral to tetragonal 
and the strong strain coupling, fpc the tetragonal direc- 
tion are the main ingredients, iia E3 Recently, it has,been 
suggested both theoreticallyo and expepmentallyEZl that 
the monoclinic phase found at the MPlS can bridge the 
tetragonal and the rhombohedral phase favoring such po- 
larization rotation. Bellaiche et al&3 used an effective 
Hamiltonian approach adapted to the alloy with param- 
eters determined from first principles calculations. The 
effective Hamiltonian used allowed strains, alloy disor- 
der and off-centering displacements to interact, but oc- 
tahedral rotations were not included. Nonetheless, ex- 
cellent agreement with experimental data was obtained 
for the temperature and composition dependence of the 
piezoelectric and structural properties in PZT, except 
that the temperature scale had to be uniformly adjusted 
downwards by approximately one third. The micro- 
scopic mechanism that makes polarization rotation oc- 
cur at modest field strengths is still somewhat uncer- 
tain though the kimpttancc of Pb chemistry has been 
widely recognized™ nr j'Bca In particular, it is not fully 
understood how the instabilities present in PbTiC>3 and 
PbZrC>3 are modified in the alloy and the extent to which 
they compete or perhaps coexist in real samples. 

Certainly, the simple supercell discussed above is a 
rather severe approximation to the disordered alloy. Even 
though Zr and Ti have the same valence, their ionic sizes 
and electronic properties differ significantly. At the very 
least, local stresses related to Zr-Ti disorder should be ex- 
pected in PZT alloys. Zr rich local regions will be under 
compressive stress so the volume available for the Pb ions 
is reduced thus lowering the local FE tendency. lonij con- 
siderations, supported by experimental evidence,Cj pre- 
dict that octahedral rotations should respond in the op- 
posite way to compression especially if the octahedra are 
stiff. Consider such rotations around [001] (C4/, sym- 
metry) and [111] (C3i symmetry). As mentioned, these 
are not true modes for the supercell, but they are zone 
boundary rotational modes that play a key role in at least 
PbZr03. Both these rotations are unstable (Fig. 2). At 
a = 7.631 ao, the energy gain for the [111] rotation is 
9 mRy (7 mRy around [001]). It is remarkable that the 
sizes of the FE and rotational instabilities are so simi- 
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FIG. 2. Variation of the energy (of 10 atom supercells) with 
"rotation" (a) around [111] and (b) around [001], for the refer- 
ence lattice parameter (circles), under compression (squares) 
and under tensile stress (diamonds) , as in Fig. jjj. 

lar for this x = 0.5 supercell, considering that they are 
also quite close, again with the FE instability lower, for 
PbZrOs though with a much higher energy scale. Pre- 
sumably, they track each other across the rhombohedral 
side of the phase diagram. This can be qualitatively un- 
derstood-iiia--terms of ionic size effects. Unlike BaTiOa or 
KNb03,Eilc3 the ferroelectric mode in these Pb based ma- 
terials is best described as a Pb off-centering with respect 
to the surrounding O with a smaller transition metal dis- 
placement in the same direction. This is clearly seen, 
e.g. in the Tis phonon eigenvector. Thus the FE insta- 
bility is controlled by the volume available to the Pb ion; 
the rotational mode also involves changes in Pb-0 bond 
length, again controlled by the same distances. 

However, a 1% compression (to a — 7.555 ao) increases 
the rotational instabilities to 11 mRy about [111] and 9 
mRy around [001]. Qualitatively, this is related to the 
stiffness of the O octahedra, especially around Ti; un- 
der pressure the octahedra size and geometry can only 
be retained by rotation, which is made at the expense of 
the relatively soft Pb-0 interaction. This driving force is 
clearly absent for the FE instability. Experimental evi- 
dence of opposite stress dependencies of rotational (R25) 
and-FE (Fi 5 ) instabilities was discussed-for PbZr03 early 
ono Stress-field response calculations^ for PbZro.sTio.5 
supercells suggest that competition between different dis- 
tortions may be important when uniaxial stress is not 
parallel to the FE distortion. Considering the differ- 
ent experimental cell volumes of the end-points (effec- 
tiva-o^PbZrOa) = 7.7883 a and a(PbTi0 3 ) = 7.5028 
ao )EJE-3, local stresses sufficient to tip the order of the 
FE and rotational instabilities seem quite reasonable in 
the disordered alloy. As noted, the FE instability oc- 
curs at the r and R points and probably over most or all 
of the zone. Thus it can occur on a very local scale in 
real space. The rotational instability is no doubt more 
localized in reciprocal space, so several octahedra must 
rotate in concert within a region, but this region may be 
realizably small, as the ZrOe octahedra may be soft as 
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in PbZrC>3 providing weak connections in the interlinked 
network. In this scenario, coexistence of rotational and 
FE distortions in the alloy is expected. The resulting 
disorder might increase response away from the MPB, 
though at the expense of maximum attainable response 
at the best composition - this would be favorable for ob- 
taining a desirable weak temperature dependence of the 
response. In any case, it is clear that rotational degrees 
of freedom are in the low energy space along with the 
FE mode, and presumably should be considered in the 
construction of effective Hamiltonians, either via explicit 
additional coordinates or rcnormalization of the existing 
FE and strain related coordinates. This could possibly 
improve the temperature scale. 
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